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Abstract 

The present study deals with dissipative future universe without big rip in context of Eckart formalism. 
The generalized chaplygin gas, characterized by equation of state p = — 4r, has been considered as a 

model for dark energy due to its dark-energy-like evolution at late time. It is demonstrated that, if the 
cosmic dark energy behaves like a fluid with equation of state p = tup; u < — 1, as well as chaplygin gas 
simultaneously then the big rip problem does not arises and the scale factor is found to be regular for all 
time. 
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Recent observations like CMB anisotropy, supernova and galaxy clustering strongly indicate that our 
universe is spatially flat and there exists an exotic cosmic fluid called dark energy with negative pressure, 
which constitutes about 70 percent of the total energy of universe. The dark energy is usually described by 
l/~, ■ an equation of state (EoS) parameter u = - , the ratio of spatially homogeneous dark-energy pressure p to 

its energy density p. A value to < — ^ is required for cosmic acceleration. The simplest explanation for dark 
energy is cosmological constant, for which u = — 1. The increasing evidence from observational data indicates 
that u) lies in a narrow strip around w = — 1 quite likely being less than this value [T] — [3] . The region where 
EoS parameter u> < —1, is typically referred to as a phantom dark energy universe. The existence of the 
region with u < — 1 opens up a number of fundamental questions. For instance, the entropy of such universe 
is negative. The dominant energy condition (DEC) for phantom fluid is violated, as a rule. The phantom 
dominated universe end up with a finite time future singularity called big rip or cosmic doomsday [4] [5] . The 
last property attracted much attention and brought the number of speculations upto the explicit calculation 
of the rest of the life-time of our universe. 



o 



X 



Soon after Caldwell [4. proposed phantom dark energy model with cosmic doomsday of future universe, 
cosmologists started making efforts to avoid this problem using u < —1 [5] — |S]. In the braneworld scenario, 
Sahni and Shtanov has obtained well-behaved expansion for the future universe without big rip problem with 
w < — 1. They have shown that acceleration is a transient phenomenon in the current universe and the future 
universe will re-enter matter dominated decelerated phase [5] . It is found that general relativity (GR) based 
phantom model encounters "sudden future singularity" leading a divergent scale factor, energy density and 
pressure at finite time t — t s . Thus the classical approach to phantom model exhibits big rip problem. For 
future singularity model, curvature invariant becomes very strong and energy density is very high near t = t s 
[ID] . So, quantum effects should be dominated for \t = t s \ < one unit of time (Early universe) [II] — [13] 
and it is shown that the an escape from the big rip is possible on making quantum corrections to the energy 
density and pressure in Fricdmann equations. 

In the framework of Robertson- Walker cosmology, Chaplygin gas (CG) is also considered as a good source 
of dark energy for having negative pressure, given as 

P=-~ (1) 

P 



where p and p are, respectively, pressure and energy density in a co moving reference frame, with p > 0; A 
is a positive constant. 

Moreover, it is only gas having super-symmetry generalization |14 | — [16 ] . Bertolami et al |17) have found 
that generalized Chaplaincy gas (CG) is better fit for latest Supernova data. In case of CG, equation ((T|) is 
modified as 

P=-Ar (2) 

pa 

where 1 < a < oo. 

For a = 1, equation ([2]) corresponds to equation ([I]). 

Other approaches have considered dissipative effects in CG models, using the framework of Elkhart theory 
[T5] . Thai et al [T5], have investigated a viscous GCG, assuming that there is a bulk viscosity in a linear 
borotropic fluid and GCG. It is found that the equation of state of GCG can cross the boundary u = — 1. 
Also in Ref. [20], it is found that a dissipative chaplygin gas can give rise to structurally stable evaluational 
scenarios. It is interesting to note that the GCG itself can behave like a fluid with viscosity in the context of 
Eckart formalism [T8]. Fabris et al [2"T] . have investigated an equivalence GCG and dust like fluid. Recently 
Cruz et al [35], have studied dissipative generalized chaplygin gas as phantom dark energy and found the 
cosmological solutions for GCG with bulk viscosity. 

The FRW metric for an homogeneous and isotropic flat universe is given by 

ds 2 = -dt 2 + a{t) 2 (dx 2 + dy 2 + dz 2 ) (3) 

where a(t) is the scale factor and t represents the cosmic time. 
The field equations in the presence bulk viscous stresses are 

— = -Up+m (5> 

a o 

where p is the effective pressure given by 

p = p- 3#£ (6) 

Here p, £ are the isotropic pressure and bulk viscous coefficient respectively. 
The energy conservation equation is given by 

p 4 + 3H(p + p)=Q (7) 

Here, and in what follows the sub indices 4 on a, p and elsewhere denote differentiation with respect to t. 
Using equations ([2]), ((4]) and (|6]), equation ([7]) leads to 

P4 + 3^( P -^-3hA=0 (8) 

a V p~ J 

In order to obtain solution of equation ([8]), we will assume that the viscosity has a power-law dependence 
upon the density 

£ = Cop" (9) 

where £o and n are constant. 

On using equation §§§ in equation ©, we obtain 

* + »*£^Lv"6 do) 

dt a dt pa 

To solve equation (fit)]) , we use the transformation p 4 = f(p) = p n+1 , accordingly equation (|10l) leads to 
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(11) 



where po = p(t ) and a = a(t ); t is the present time. 

In the present model, it is assumed that the dark energy behaves like GCG, obeying equation $Z§ as well as 

fluid with equation of state 

p = Ljp (12) 

with ui < — 1 simultaneously. 

From equations ([2]) and (12), we obtain 
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So, evolution of equation ([TB")) at t = to leads to 



A = -uj p ' 



(13) 



(14) 



with uiq = u>(to). 

Using equation (fH)) . equation ([TTj) leads to 



P = Pa 



a \ =a-3«o) 
- Wo + (l+u;o)(^ 



In the homogeneous model of universe, a scalar field 4>{t) with potential V(4>) has energy density 



and pressure 



P^ = 2^4 + V{4>) 



P4> = \<& ~ V{cf>) 



K = p<p + p<t> 



Equation ([16]) and dTTJ) lead to 

Using equations @, (12) and (fTl)) , equation (fig)) reduces to 
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Equation ([15]) and ([T<J]) lead to 
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From equation (|20|) . it is clear that for 1 + ujq > 0, 4>\ > 0, giving positive kinetic energy and for 1 + loq < 0, 
04 < 0, giving negative kinetic energy. 1 + ujq > and 1 + wo < are representing the case of quintessence 
and phantom fluid dominated universe respectively. Similar result are obtained by Hoyle and Narlikar in C 
-field with negative kinetic energy for steady state theory of universe [23] . 
Now, from equation ^ and ([15)) , we obtain 
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(21) 



where \ujq\ 



-td, Hq = 100hkm/s Mpc, present value of the Hubble's parameter and S7q 
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Equation (|2"Tj) may be written as 
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(22) 
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Figure 1: The plot of scale factor (a) and horizon distance (djj) versus time (t). 
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Neglecting the higher powers of ( ij, | ( ^fy ) " , equation (1221) leads to 
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= V^o-ffo|wo| 2(1+a) 



a{\ - |wo|) ( ao \ Q ' <1 - +3 ^ ) 
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Integrating equation (|23[) . we obtain 



a(t) = 
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(23) 
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From equation (|24p . it is clear that as t —¥ oo, a(t) — > oo which is supported by recent observation of Supernova 
la [231 US] and WMAP [25] [57] . Therefore the present model is free from finite time future singularity. 
Now the horizon distance is obtained as 



d H {t) = 



3(l + o)o(t) ( 2(l + a)|wo| \ w' 



a(l-3£ )a V" + (a + 2)|a7 | 
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(25) 

Fig. 1 depicts the variation of scale factor (a) and horizon distance (djj) versus cosmic time, as a represen- 
tative case with appropriate choice of constants and other physical parameters. From equation (|2"4"|) and ([23]) 
it is clear that du(t) > a(t) i. e. horizon grows more rapidly than scale factor which is clearly shown in Fig. 



In this case, the Hubble distance is given by 
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(26) 
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Figure 2: The plot of Hubble distance (H 1 ) versus time (£). 



Equation (|26|) is showing the growth of Hubble distance (H 1 ) with time such that H 
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as t — !• oo. This behaviour of H~ l is clearly depicted in Fig. 2. Thus in present case, the galaxies will 
not disappear when t — > oo, avoiding big rip singularity. Therefore, one can conclude that if phantom fluid 
behaves like GCG and fluid with p = ujp simultaneously then the future accelerated expansion of universe 
will free from catastropic situation like big rip. Equation (|15|) may be written as 



P = Po 



|wo| + (1 - |wo|) 






*(l-3« ) 



(27) 



From equation (|27l) . it is clear that as t — > oo, p — > pq\ll> ( )\ t +^ > p (since t —} oo, a(t) —} oo). Thus one can 
conclude that energy density increases with time, contrary to other phantom models having future singularity 
at t = i s [4j [IT] . Based on la Supernova data, singh et al [28] have estimated ujq for model in the range 
—2.4 < loq < —1.74 up to 95 percent confidence level. Taking this estimate as an example, with a = 1, Poo 
is found in the range 1.31po < Poo < l-54po and with a = 2, p^ is found in the range 1.44po < poo < 1.78po 
and so on. This does not yields much increase in energy density as t — > oo but if the future experiments 
supports large value of |wo| then p^ will be high. 

It is interesting to see here that present model is derived by Bulk viscosity in the context of Eckart 
formalism 18 . We see that big rip problem does not arises in the present model. In Refs [TD] — [13]' f° r 
models with future singularity escape from big rip, is demonstrated using quantum corrections in the field 
equations near t — t s . However the present study deals phantom cosmology with accelerated expansion 
without catastrophic situations using classical approach. It is also seen that the model of future universe 
presented by Srivastava [29] is a particular case of the model presented in this letter. 
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